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We present a very high statistics study of D and D s semileptonic 
+^ decay form factors on the lattice. We work with MILC Nf = 2 + 1 

T< lattices and use the Highly Improved Staggered Quark action (HISQ) 

Ph for both the charm and the strange and light valence quarks. We use 

r-| both scalar and vector currents to determine the form factors fo(q 2 ) and 

f+{q 2 ) for a range of D and D s semileptonic decays, including D — > tt£v 
t-h and D — > Ktv. By using a phased boundary condition we are able to 

^ tune accurately to q 2 = and explore the whole q 2 range allowed by 

■^j- kinematics. We can thus compare the shape in q 2 to that from experiment 

and extract the CKM matrix element \V r A. We show that the form factors 
are insensitive to the spectator quark: D — > Ktv and D s — > i] s iv form 
q factors are essentially the same, which is also true for D — > iriu and 

CN D s Kiv within 5%. This has important implications when considering 

the corresponding B/B s processes. 
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1 Introduction 



Lattice QCD is an excellent tool for calculating strong interaction effects from first 
principles and can provide accurate phenomenology not available with any other 
method. Particularly important and interesting applications of this method are 
the determinations of various heavy meson weak decay matrix elements that are 
key to constraining the vertex of the Unitarity Triangle derived from the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix. This provides a stringent test of the self-consis- 
tency of the Standard Model. 

However, calculations of these matrix elements must not be seen in isolation. This 
is only one part of the HPQCD program to calculate meson spectra, decay constants 
and other QCD observables fully non-perturbatively (see e.g. the calculation of fjj g 
in pQ or J/ip mass, leptonic width and radiative decay rate to rj c in [2]). The Highly 
Improved Staggered Quark (HISQ) formalism enables us to keep the discretization 
errors small and to treat charm quarks the same way as light and strange quarks, 
which reduces the systematic errors. 

A simple recipe for doing a Lattice QCD calculation is: 

1. Generate sets of gluon fields from Monte Carlo integration of the QCD path 
integral (including effects of u, d and s sea quarks) 

2. Calculate averaged hadron correlators from valence quark propagators 

3. Fit the correlators as a function of time to obtain masses and matrix elements 

4. Determine a and fix m q to get results in physical units 

5. Extrapolate to a = and physical light quark mass for real world 

However, in this paper we don't go into details of Lattice QCD and concentrate on 
the results. The paper is organized as follows: In Section [2] we briefly explain how we 
calculate the form factors on the lattice and fit the results. We give our results for the 
form factors in Section [3] and talk about the ^-expansion and continuum and chiral 
extrapolations in more detail in Section |4| In Section [5] we explain how we extract 
V cs and give our result and summarise in Section [6j 



2 Form factors on the lattice 

Semileptonic form factors are 3-point amplitudes in Lattice QCD. To get information 
of both form factors, /o(<? 2 ) an d f + (q 2 ), we calculate scalar and vector currents and 
the corresponding 2-point correlators for the mesons - see Fig. [T] The scalar current 
is 

{K\S\D) = fr K {q 2 ) Ml ~ M2K , (1) 

m 0c - m 0s 
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Figure 1: 2-point and 3-point correlators. In the 2-point correlator a kaon is created 
at time t' and annihilated at time 0. In the 3-point correlator the two mesons are 
time T apart and a scalar or vector current J is inserted at time t (0 < t < T). One 
of the quarks can be given a momentum p to explore the full q 2 range. 

where q^ = Pd ~ Pk (using D — > Klv as an example), and the vector current can be 
written as 



(K\W\D) = f^ K (q 2 " 



Pd+Pk o ? 

q z 



+ fo D ^ K (Q 2 ) ° a V - (2) 



<7 2 



Note that this guarantees that /o(0) = /+(0). See also our earlier calculation of 
semileptonic form factors in [3]. We use 3 ensembles of MILC asqtad Nf = 2 + 1 
lattice configurations (sets 1, 2 and 4 in [2]). 

From experiments we get the differential decay rates, e.g. 



= Gjrpj; 

dq 2 24vr 3 



V cs \Vr K (lT (3) 



for D — > Klv. From this one can determine \V CS ■ f+^ K {q 2 )\, but one needs either 
f+{q 2 ) from theory or V cs from unitarity to determine the other. 



2.1 Fitting the Lattice QCD results 

We fit the lattice 2-point and 3-point correlators simultaneously, as a function of 
time i! and meson separation T (see Fig. [TJ to estimate correlations between all fit 
parameters. We use multi-exponential fits (up to 5 exponentials) to reduce systematic 
errors from the excited states. The fit parameters are constrained by Bayesian priors. 
More details about the type of fits that we use can be found in [3]. 

We have very high statistics, of the order of 100000 correlators, which gives us 
very small statistical errors. By giving one of the quarks a momentum p using twisted 
boundary conditions [1] , as shown in Fig. [TJ we can tune accurately to q 2 = or choose 
any q 2 value in the allowed kinematical region. This enables us to study the whole 
physical q 2 range. 
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3 Semileptonic form factors 



We have calculated the form factors fo and /+, as a function of q 2 , for various D and 
D s semileptonic decays. Let us divide the results into two groups: charm to strange 
decays and charm to light decays. 

3.1 Charm to strange decay: D — > K£v, D s — > rj s iu and D s — > 



s I c \'")+ \ 



K ^ \fl 




spectator quark 



Figure 2: On the left: D — > Ktv and D s — > rjsiu decays are the same except for the 
spectator quark; On the right: Kinematics of the D s — > cf)lv decay. 

D — > Kiv and D s — > i] s £u both have the same charm-strange current, as they 
are both charm to strange decays. rj s is the pseudoscalar ss meson, so both decays 
are pseudoscalar to pseudoscalar decays. Note that i] s is not a physical meson, but 
can be easily calculated on a lattice. The difference between these two decays is the 
spectator quark, light vs. strange, as illustrated in Fig. [2j 

Our results on different lattice ensembles are shown in Fig. [3j along with the 
final result after the continuum and chiral extrapolation. The q 2 dependence is well 
understood qualitatively - / + rises more steeply than f , as /+ is governed by the 
vector meson M^. in the pole mass parameterization, whereas /o is governed by the 
scalar meson Md* . The continuum and chiral extrapolation will be discussed in more 
detail in Section |4j The differences between the coarse and the fine lattice results are 
small, i.e. the discretisation effects are very well under control. 

Note that the shapes of the form factors do not depend on the spectator quark: the 
form factors for these two decays, D — > Klv and D s — > r] s £i>, are the same within 3% 
and even closer when one moves away from q 2 = 0. Comparing the decay constants 
of the two mesons, fo and fr> s , one might expect a change of about 15% when going 
from a light spectator quark to a strange quark. However, this does not appear to be 
the case for form factors. 

On the lattice one can also calculate form factors for a pseudoscalar to a vector 
meson semileptonic decay. D s — > cf)iv is a charm to strange decay like D — > Klv and 
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Figure 3: D — > Klv and D s — > r\ s tv form factors. The form factors are very insensitive 
to the spectator quark. 

D s — > 7] s [v. As is a vector meson, there are now more form factors than in the 
pseudoscalar to pseudoscalar decay. The general form for a pseudoscalar to vector 
matrix element is 

W, e)|V - A"|£>.(p)> = A/ f £ T" t u PlAv{q 2 ) - (M Ds + M^Atf) 

+ M € '1 M (P + P')^(g 2 ) + 2M^gM 3 (g 2 ) - 2M^gM (g 2 ), (4) 

where 

Mi 2 ) = "^^^ - ^%^W> md *«» = -V *- (») 

Here p is the momentum of the D s , p' is the momentum of the <fi and e is the polari- 
sation vector. V(q 2 ) is the vector form factor and A (q 2 ), A 1 (q 2 ), A 2 (q 2 ), A 3 (q 2 ) are 
axial vector form factors. Note that only three of the axial vector form factors are in- 
dependent. We can extract the different form factors by choosing the right kinematics 
- see [5] for more details. The diagram on the right in Fig. [2] shows the kinematics 
of this decay: 9 e is the angle between the e + and the W and 9k is the angle between 
the kaon and the 0. 8 e is measured in the rest frame of the W and 9k is measured in 
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Figure 4: D s — > cf)iv form factors, defined in [3]. 



the rest frame of the 0. Our results, the form factors and angular distributions, are 
plotted in Figs. [4] and [5j The 'boxes' in Fig. [5] show the experimental results from 
BaBar [8] for the angular distributions - we see good agreement between theory and 
experiment. 



3.2 Charm to light decay: D — > irii/ and D s — > Klv 

We calculate two charm to light semileptonic decays, D — > iriu and D s — > Kiv. In this 
case, both of the decays are experimentally accessible. Our results on different lattice 
ensembles, along with the final result after the continuum and chiral extrapolation 
(see Section [4]), are shown in Fig. [6] The conclusions are very similar to the charm 
to strange decay: Again, the dependence of the form factors on the spectator quark 
mass is very mild - going from a strange spectator quark to a light quark changes the 
form factors by less than 5%. 
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Figure 5: D s — > <\)lv decay rates in q 2 bins and angular distributions. The angles are 
defined in Fig. [2] 

4 The z-expansion 

For continuum and chiral extrapolations we use the z-expansion |6|: First we remove 
the poles (using D — >• K decay here as an example) 



i 



ML 

U s0 ' 



M 2 



and convert to z space 



t + -q* 



t+ = (M D + M K ) 



f+(q 2 ) (6) 



(7) 



t+ - q 2 + 

The transformation of the complex q 2 plane to the z plane is sketched in Fig. [7j This 
is useful, as the form factors / in the physical region can be described by a simple 
power series in z: 

/^( 2 ) = ^„(a)/, f^ K (z) = cn(a)z n , c = b . (8) 

n>0 n>0 

We let the fit parameters b n , c n depend on lattice spacing and sea quark masses. In 
the end we take a = and m q = to get the result in the physical limit. Fig. [7] 

shows one of the fits, for D — > Klv. 
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Figure 6: D — y Triu and D s — y Klv form factors. The form factors are not sensitive 
to the spectator quark. 

5 Extracting V cs 

One of the goals is to determine the CKM matrix elements and over-constrain the 
unitarity triangle by providing accurate input from theory. The lattice calculation 
does not know about V cs , but experiments do. By combining the lattice calculation 
of the D — y Klv form factor f+{q 2 ) with experimental results we can extract V cs 
without the need to assume unitarity of the CKM matrix. By integrating the lattice 
form factors over the experimental q 2 bins, i.e. integrating Eq. [3] to get the rate for a 
given bin and then taking the experiment to lattice ratio gives us V 2 S for that given 
bin. Our results are shown in Fig. [8j We do this using CLEO [7j, BaBar [8], Belle [9] 
and BESIII (preliminary, [10]) results and fit a constant to these V 2 S values, including 
bin to bin correlations from lattice calculations and experiments in the fit. Including 
different sets of experimental results does not change our value for V cs , as can be seen 
in Fig. [9j Our best, preliminary value is V cs = 0.965(14). This is consistent with 
V cs = 0.97344(16) from PDG [12] . that is calcuated by assuming unitarity. Using our 



best value for V cs we plot the D — y Klv decay rates in q 2 bins in Fig. 10, which shows 
excellent agreement between theory and experiments. D s — y <f)£i> decay rates, plotted 
in Fig. [5j also show excellent agreement. 
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Figure 7: On the left: D — > Ktv form factors in z space and the fit. On the right: 
Transformation of the complex q 2 plane to the z plane. 



6 Summary 

This very high precision Lattice QCD calculation provides D and D s meson semilep- 
tonic decay form factors from first principles. We study the full q 2 range and different 
semileptonic decays (different daughter mesons). We determine the form factors to 
better than 3% accuracy (better than 2% in the case of D — > Ktv). The D/D s form 
factors are very insensitive to the spectator quark and this is expected to be true for 
B/B s as well. We calculate decay rates in q 2 bins to compare with experiments and 
extract V cs . We get V cs = 0.965(14) (preliminary), which gives very good agreement 
with experiments, i.e. the shape of the form factor f+~^ K calculated in Lattice QCD 
agrees with experimental results. 
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Figure 8: Ratio of experimental to lattice results in each q 2 bin for D — > K£v, i.e. 
\V CS \ 2 extracted from that bin directly. The experimental results are from [71 181 191 [TO] . 
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Figure 9: Our determination of V cs from this Lattice QCD study combined with 
various sets of experimental results. D — > Kiv experimental results are from [3 El 13 
[TO] . The decay constant fo a is from [lj and leptonic decay rates are from [TT] . 



16 r p$ 



'a 

A 



Oh 



18 
16 
14 
12 
10 
8 
6 
4 
2 
0, 



i i i i i r 



T 



O Lattice 

□ CLEO 

<> BaBar 

A Belle 

> BESIII prelim. 

V from combined fit 



J I I I I I I L 



0.2 0.4 0.6 o.: 



1.0 1.2 1.4 1.6 2 Total/10 

" max 

2 2 4 

q [GeV /c ] 



Figure 10: Decay rates in q 2 bins for D — > Kiv. The lattice numbers include our 
result for V cs . The experimental results are from El El ITU] . 
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